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[57] ABSTRACT 

Phase of each transducer element of a high intensity focused 
(HIFU) transducer array is controlled to compensate for 
phase change introduced by varying velocity through dif- 
fering tissue along a path towards a treatment volume. The 
echo of a specific harmonic of ultrasound pulses of moderate 
intensity and less than HIFU intensity are used to measure 
the propagation path transit time of each HIFU transducer 
element that will converge in a treatment volume through 
nonhomogeneous tissue. The moderate intensity is outside 
the linear region between molecular velocity fluctuations 
and pressure fluctuations. Thus, specific harmonic echoes 
are distributed in all directions from the treatment volume. 
Temporal delay in the specific harmonic echoes provide a 
measure of the propagation path transit time to transmit a 
pulse that will converge on the treatment volume. 

25 Claims, 10 Drawing Sheets 
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METHOD FOR DETERMINING PHASE 
ADVANCEMENT OF TRANSDUCER 
ELEMENTS IN HIGH INTENSITY FOCUSED 
ULTRASOUND 

FEDERAL FUNDING STATEMENT 

This invention was made with government support under 
grant number N0001496- 1-0630 awarded by the United 
States Navy. The government has certain rights in the 
invention. 

BACKGROUND OF THE INVENTION 

This invention relates generally to the use of ultrasonics 
in medical technology applications, and more particularly to 
a method and apparatus for applying high intensity focused 
ultrasound through nonhomogeneous volumes. 

Sound waves that have a frequency greater than approxi- 
mately 20 kHz are referred to as ultrasound. Much of the 
research and development in the field of medical ultrasound 
relates to diagnostic applications and therapeutic applica- 
tions. Medical diagnostic ultrasound systems are useful for 
generating images of anatomical structures within a patient's 
body. The images are obtained by scanning a target area with 
waves of ultrasound pulses. In therapeutic ultrasound 
applications, high intensity ultrasound pulses are transmitted 
into a target area to induce changes in state of the target. 
High intensity ultrasound pulses induce changes in state 
through thermal effects (e.g., induced hyperthermia) and 
mechanical effects (e.g., induced cavitation). 

In medical ultrasound applications, ultrasound pulses are 
delivered to a patient by an ultrasound transducer. To obtain 
the ultrasound pulses, electronic signals are input to the 
transducer. The transducer then converts the electrical sig- 
nals into ultrasound pulses which arc transmitted into the 
patient's body as an ultrasound beam. Such ultrasound beam 
is absorbed, dispersed, and reflected. Diagnostic applica- 
tions take advantage of the reflected ultrasound energy, 
which is analyzed and processed to generate image and flow 
information. Therapeutic applications take advantage of the 
absorbed ultrasound energy to change the state of a target 
area. Dispersion is an aspect of ultrasound technology that 
poses challenges to be overcome. Specifically in therapeutic 
applications, it is desired that ultrasound energy cause a 
change of state at the target area and not adversely impact 
other volumes within the patient. Refractive dispersion 
diminishes the therapeutic effect of ultrasound energy and 
causes the ultrasound energy to be absorbed at unintended 
areas. Accordingly there is a need for accurately focussing 
the HIFU beam. 

In a hyperthermia HIFU treatment the ultrasound beam is 
highly focused. The beam intensities increase along the 
beam path from the transducer to the target area. At the 
target area very high temperatures can be induced. The 
absorption of the ultrasound energy at the target area induces 
a sudden temperature rise (e.g., tens of degrees centigrade 
per second). This temperature rise causes coagulation or 
ablation of target area cells. Accordingly, hyperthermia 
HIFU treatments can cause damage to an internal lesion. It 
is desirable that such damage occur without harming inter- 
mediary cells between the transducer and the target area. 

HIFU treatments are desired for rapid heating of human 
tissues to arrest bleeding and to ablate tumors. Typically, a 
target area is located deep in the body and the ultrasound 
energy is delivered through the skin. Using an intensity in 
excess of 1000 W/cm 2 spatial peak continuous wave 
(SPCW), the target area is heated at a rate of approximately 
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25° C./second. To limit the rapid heating effects to the 
desired target area, and spare the intervening tissue, a large 
aperture transducer is used. Along the path to the target area, 
however, the ultrasound energy travels through fat and 

5 muscle, (i.e., the ultrasound travels through a nonhomoge- 
neous medium). The ultrasound travels at different speeds 
through the different materials resulting in refractive disper- 
sion of the ultrasound pulses. This causes difficulty in 
focusing to the desired target area. This invention is directed 

10 to improving focus of high intensity ultrasound through 
nonhomogeneous tissue so that (i) thermal effects in inter- 
mediary tissue are controlled and harmless, and (ii) an 
elleclive therapeutic dose is delivered to the target volume. 
The intensities described in this application are specified 

35 as spatial peak continuous wave (SPCW), spatial peak 
temporal average (SPTA), or spatial peak temporal peak 
(SPTP). Referring to FIG. 1A, a SPCW waveform 21 is 
shown. The dotted portion 23 indicates negative pressures 
which arc cut off. Conventional diagnostic ultrasound appli- 

20 cations use ultrasound pulses with intensities of up to 100 
mW/cm 2 spatial peak temporal average (SPTA) as exempli- 
fied by waveform 25 of FIG. IB. The dotted portion 27 
indicates negative pressures which are cut off. Referring to 
FIG. 1C, a SPTA waveform 29 is shown. Each of these 

25 waveforms 21, 25, 29 are for a 1 MHz ultrasound waveform. 

SUMMARY OF THE INVENTION 

According to the invention, the phase of each transducer 

3(J element of a high intensity focused ultrasound (HIFU) 
transducer array is controlled to compensate for phase 
change introduced by differing tissue along a path towards 
a treatment volume. 

Conventional diagnostic ultrasound applications use 

35 ultrasound pulses with intensities of up to 100 mW/cm 2 
spatial peak temporal average (SPTA). HIFU applications 
use ultrasound pulses with intensities of about 1000 W/cm 2 
SPCW (i.e., about 10,000 times greater). According to an 
aspect of this invention, an aiming aperture array momen- 

40 tarily generating 10,000 mW/cm 2 spatial peak temporal 
peak (SPTP) at the target volume (i.e., about the maximum 
diagnostic intensity and about 1000 times less than HIFU 
intensities) is used for determining path length for each 
element of a transducer array between the transducer head 

45 and treatment volume. In diagnostic applications using 
intensities of 100 mW/cm 2 SPTP sinusoidal pressure fluc- 
tuations are approximately 0.6 atmospheres. Thus, the high- 
est instantaneous pressure in tissue is between 0.4 atm and 
1.6 atm. At these intensities, a linear relationship exists 

50 between molecular velocity fluctuations and pressure fluc- 
tuation. At 300 mW/cm 2 SPTP the pressure variation is 1.0 
atm. Thus, the highest instantaneous pressure at such inten- 
sity is between 0.0 atm and 2.0 atm. Within this pressure 
range, there still is a linear relationship between molecular 

55 velocity fluctuations and pressure fluctuations. The 300 
mW/cm 2 SPTP intensity is called the threshold intensity. 
Further increases in intensity cannot drive the lower limit 
negative. As a result, at higher intensities the pressure wave 
becomes clipped, introducing harmonics into the pressure 

60 fluctuations. Echoes from the lower intensity regions will 
contain only the fundamental frequency, whereas echoes 
from the target volume (i.e., a high intensity region) will 
contain the fundamental frequency and its harmonics. 
Because the ultrasound transducer is sensitive to the odd 

65 harmonics, a specific harmonic, such as the third harmonic- 
echo of an ultrasound pulse, is used to measure propagation 
path transit time. A moderate power intensity, greater than 
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the threshold intensity (300 mW/cm 2 SPTP) and less than Application of the ultrasound energy causes molecular 
HIFU intensity (1000 W/cm 2 SPCW), is used. The propa- velocity fluctuation and pressure fluctuation along a path 
gation path transit time is measured for each transducer from a source of ultrasound energy to the treatment volume, 
element's pulse that will converge in a treatment volume There is a threshold intensity level of the applied ultrasound 
through nonhomogcneous tissue. The moderate power level 5 energy, above which a relationship between the molecular 
10 W/cm 2 SPTP is not enough to be considered HIFU and velocity fluctuations and the pressure fluctuations becomes 
will be affected by refractive aberration by, for example, nonlinear. The method of treatment includes, at one step, 
subcutaneous fat, along the path between the transducer and identifying a treatment volume to be treated. At another step 
the treatment volume. Although such aberration will defocus configuring a firsl lrans ducer array to define a first aperture, 
the beam, decreasing the intensity at the treatment volume, ]Q ^ first transducer ^ formed 5 a luralit of trans . 
the moderate intensity must be above the threshold value so , . . J a j • a ^ 
that harmonics are generated in and radiated from the du T^Tt ^ . ^ " Ta y» confl »l^ 
treatment volume. Because the 100 W/cm 2 SPIT power <° def \ ne , the first a ? erUire ' ^trasound energy is transmitted 
level is above the linear region between molecular velocity from . lhe f,rst transducer array at a prescribed frequency and 
fluctuations and pressure fluctuations, the echoes at the at a first intensity greater than the threshold intensity. Echoes 
specific harmonic used are distributed in all directions from 15 of the transmitted ultrasound energy are received at the first 
the treatment volume. These harmonics are coherent with transducer array and at a second transducer array having a 
the fundamental frequency. Thus, the temporal delay in the plurality of transducer elements. The echoes include ultra- 
specific harmonic echoes observed by each transducer sound energy at the prescribed frequency and at harmonics 
element, such as the third harmonic echoes, provides a of the prescribed frequency. The phase of the echoes are 
measure of the propagation path transit time that will be 20 compared at selected harmonics, such as the third harmonic, 
experienced by the corresponding portion of the transmitted for each element of the first and second transducer arrays. 
HIFU ultrasound, if it is to converge on the treatment Steering and/or focus delays for one or more elements of the 
volume. transducer arrays then is adjusted to alter the phase, and 

According to a method of this invention, a diagnostic optionally the amplitude, of the received echoes. Steps of 

ultrasound imaging transducer array is used to identify and 2 5 transmitting, receiving, comparing and adjusting are 

locate a treatment volume to be treated with HIFU energy. repeated until each element of the first and second trans- 

A high power test aperture including imaging transducer ducer array is properly aimed at the treatment volume. Such 

elements then is used to focus on the treatment volume and adjusting step serves for a phase inversion process in which 

deliver a super threshold test pulse. At another step echoes respe ctive propagation transit times are applied as an 

at a specific harmonic, such as the third harmonic, are 3Q advance time t0 each corresponding transducer element, 

received by all elements of an enlarged high intensity While phasc inversioD fa appliedj ultrasoun d energy is 

aperture. Phase is measured for each clement. A phase lransmitted from the first an d second transducer array to the 

inversion delay then is computed for each element. At treatment volume at a second intensity for achieving desired 

another step the enlarged high intensity aperture is used to ultrasonic medical therapy of tne treatment volume. The 

apply high intensity ultrasound energy focused to the treat- 35 second intensity at the treatment volume is greater than the 

ment volume. ^ rst intensity and is a HIFU intensity. In some embodiments 

According to another aspect of the invention, pitch tne first and second transducer arrays are portions of one 

between transducer elements is prescribed to be finer than transducer array. In some embodiments the steering and 

the spatial frequency of undulations of subcutaneous fat. As focus delay adjustment only is performed for elements of the 

ultrasound energy passes through subcutaneous fat refrac- 40 ^ond transducer array 

live aberration occurs Such subcutaneous fat typically is " Accofdin t0 one advaDtage of the invention, HIFU 
undulated in a wave-like shape bavmg different thicknesses ^ are lransmjUed tlm h non homogenec)us tissue with- 
al different areas. By spacing the transducer elements at a ou( advcrse| heatj tjssue in lhe atn belween lne ultra . 
pitch finer than the spatial frequency of these undulations lransducer and the lreatment volume . and other 
the vanations in transit time along each of the ultrasound 45 a an£j advant of the invention will be better 
trajectory paths from each individual transducer element to underslood b refere „ ce lo the following detailed descrip- 
the commmon treatment volume can be compensated by (i()n [aken fa conjunction with the accompa n y ing drawings, 
electronic delay lines, one located between each transducer 

element and the common HIFU transmitter, so that the BRIEF DESCRIPTION OF THE DRAWINGS 

transit time along each path consisting of the delay line, the 50 F]GS 1A _ C ^ waveforms of a L MHz ultrasound signal 

transducer element and the Ussues (interposed fat, muscle dm inlcasilies exem plifying SPCW amplitude, SPTP 

and other tissues), between the common HIFU transmitter ... , to , CDTA ort ^rf»^. 

" . ... I amplitude and SPTA amplitude; 
and the common treatment volume is identical. 

. t c t . . . . FIG. 2 is a block diagram of an ultrasound system 

According to another aspect ot the invention, the power ... 6 ... . . ' 

... - , a ' e tU . . 1 according to an embodiment of this invention; 

delivered to select elements of the transducer array is 55 & 

adjusted to improve the convergence of power on the FIG. 3 is a diagram ultrasound parameter control hnes for 

treatment volume. The select elements arc the transducer the transducer array of FIG. 2; 

elements terminating a path from the treatment volume FIG. 4 is a diagram of an ultrasound transducer array 

where the ultrasound has been subject to excessive attenu- layout for the transducer array of FIG. 2; 

ation. 60 FIGS. 5A-B are side view diagrams of the ultrasound 

According to another aspect of the invention, select transducer array receiving echo wavef rants from a treatment 

elements of a transducer array overlaying a patient's ribs are volume; 

deactivated to avoid bone healing, without degrading the FIG. 6 is a first schematic view of a HIFU treatment to an 

effectiveness of other transducer elements emitting ultra- organ; 

sound toward the treatment volume. 65 FIG. 7 is a second schematic view of the HIFU treatment 

During a HIFU treatment, ultrasound energy is applied of FIG. 5 showing ultrasound energy applied between the 

through a non homogeneous mass to a treatment volume. intercostal openings of a patienl\s. rib cage; 



01/20/2004, EAST version: 1.4.1 



6,042,556 

5 6 

FIG. 8 is a How chart of a method for focussing the by the transmitter. The focus control defines respective 

transducer array of FIG. 4 according to an embodiment of delays for the firing of the active transducer elements. The 

this invention; steering control defines a beam tilt resulting from timing 

FIGS. 9A-B are more detailed flow charts of a method for delays. The aperture may vary for transmission and reccp- 

focussing the transducer array of FIG. 4 according to an 5 lion. The focus and steering preferably are the same for 

embodiment of this invention; transmission and reception. Preferably, the receiver aperture 

FIG. 10 is a diagram of echo signal frequency magnitudes includes a11 transducer elements with the data for each 

over time for a transmit signal having a prescribed frequency channel processed separately. 

(three cases are shown for three different intensities of the Durin S the calibration process of this invention, the 

transmit signal)- 10 P rocessin g anc * control unit 30 (see FIG. 3) adjusts the 

. c ... , . . controls automatically to accurately focus the ultrasound 

FIG. 11 is a diagram of third harmonic frequency mag- . • ' ' 1jC 

. - ? . . .. . I j burst on the desired treatment volume 16. Pulse echoes of 

nitudes over time for the same transmit signals; and . . . n . . . r . . P 

„ „ . - . - , . r..^^ <rt ultrasound bursts are reflected back from various tissue of 

FIG. 12 is a diagram of the ratio of the curves in FIGS. 10 ^ paljem tQ lhe u|lrasound lransUucer array 12 . ^ 

an is received ultrasound energy is converted into electronic 

DESCRIPTION OF SPECIFIC EMBODIMENTS signals. Aperture, focus and steering parameters determine 

Overview the timing and adjustment of the electronic signals. The 

FIG. 2 shows a block diagram of an ultrasound system 10 electronic signals arc sampled at the receivers 28. The 

for performing high intensity focused ultrasound ('HIFU*) received signals arc processed by the processing and control 

according to an embodiment of this invention. A HIFU beam 20 unit 30. The processing unit 30 performs signal processing 

14 is emitted from a transducer array 12 and passes through analysis. In addition the processing and control unit 30 

an acoustic coupling pad 13 into a patient P to therapeuti- embodies conventional image processing and video process- 

cally act upon the patient. The ultrasound energy is in part ing methods for generating images to be displayed on a 

absorbed by tissue or other physical matter along the path of display device. Furthermore, the processing and control unit 

the beam 14. FIG. 2 shows various wave fronts 17 of the 25 30 implements software embodying a method for effectively 

ultrasound beam 14. The absorbed energy at the focus of the focusing the transmitted ultrasound burst according to an 

ultrasound beam causes thermal and mechanical state embodiment of this invention. Specifically, the processing 

changes of the tissue and other physical matter in the focus and control unit 30 provides phase error inversion adjust- 

region of the patient. In particular, the HIFU beam 14 ments to the delay lines 24. The processing unit 30 is formed 

induces hyperthermia and/or tissue necrosis at a desired 3U by one or more microprocessors, digital signal processors, 

treatment volume 16. Very high temperatures (e.g., 80° C.) multiprocessors or other known processing unit structure, 

are quickly induced (e.g., tens of degrees centigrade per The program implementing the method of this invention 

second increases) at the treatment volume causing ablation adds novel configurations to the processing unit 30. 

of treatment volume cells or necrotization of an internal FIG. 4 shows a transducer array 12 layout according to an 

lesion. The ultrasound energy may be dispersed by refrac- 35 embodiment of this invention. The transducer array includes 

tion along its path to the treatment volume 16. This invention an imaging array portion 36 and a HIFU array portion 38 

is directed toward accounting for such refractive dispersion concentrically surrounding the imaging array portion 36. 

and compensating for the dispersive effect of refraction to The imaging array portion 36 includes a plurality of trans- 

effectively focus the HIFU beam 14 upon the treatment ducer elements 34a through 34£. In one embodiment the 

volume 16 and achieve desired therapy. 40 imaging array portion 36 is formed by a linear transducer 

The ultrasound system 10 includes a user interface 18, a array. In other embodiments the imaging array portion 36 is 
system controller 22, the transducer array 12, adjustable formed by a circular, triangular or multidimensional trans- 
delay lines 24, transmitters 26, receivers 28 and a processing ducer array (e.g., 1.5 dimensional array or a 2 dimensional 
and control unit 30. In one embodiment there is a transmit array). The HIFU array portion 38 includes a plurality of 
channel and a receive channel, or a transmitter 26 and a 45 transducer elements 32a~32n. In one embodiment the HIFU 
receiver 28 for each element 32, 34 of the transducer array array elements 32 and the imaging array elements 34 are of 
12. The system controller 22 serves to control the timing of the same size. In another embodiment, as shown in FIG. 4, 
the system and triggers the transmitters 26 and receivers 28. the HIFU transducer elements 32 have a larger surface area 
For the transmit function, the system controller generates than the imaging transducer elements 34. In an alternative 
electronic signals at desired frequencies which are input to 50 embodiment the imaging transducer array portion 36 and 
the transmitters 26 to trigger an ultrasound burst. For a HIFU transducer array portion 38 are separate transducer 
receive function the system controller triggers the receivers arrays having known locations relative to each other. 
28 to sample the ultrasound echoes. FIGS. 5A and 5B show the transducer array 12 along with 

The transmitters 26 drive the respective elements of the spherical wavefronts of the ultrasound echoes received from 

transducer array 12. The receivers 28 sample the returning 55 the treatment volume 16. Note that the larger transducer 

echoes. The transmitted ultrasound burst is controlled using elements 32 are tilted toward the treatment volume. This is 

delay, focus and steering controls. The received echoes are preferable so that the spherical wavefronts 101ff-101d cross 

controlled using delay, focus and steering controls. In one a given transducer element 32 at a given angle. The tilt may 

embodiment the transmitter and receiver share a common vary for each of one or more elements 32. Phase adjustments 

steering control. The controls are implemented by the pro- 60 then are made to compensate for refractive aberration, 

cessingand control unit 30 and the adjustable delay lines 24. Method for Focussing High Intensity Ultrasound Beam 

FIG. 3 shows the delay lines 24 and the transducer array 12. High intensity ultrasound is used to oblate tissue, to 

Pulses from the transmitters 26 travel along a respective coagulate blood in tissue, and to kill micro-organisms and 

delay line to a corresponding transducer element 32, 34. other infective agents. The volumes treated with the high 

Aperture is a control of which respective transducer ele- 65 intensity ultrasound typically are deep within normal organs 

ments 32, 34 of array 12 are active. Thus, the aperture and tissues. A high intensity focused ultrasound beam deliv- 

control determines which transducer elements are activated ers high power to a small treatment volume (e.g, cubic 
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millimeter volume range), while sparing treatment from the 
surrounding tissue. One measure of treatment is referred to 
as exposure. Exposure equals the intensity of the ultrasound 
beam times the absorption rate. Exposure is to be maximized 
in the treatment volume and minimized in all surrounding 
tissues. The surrounding tissues at greatest risk from treat- 
ment harm are those superficial tissues that are between the 
treatment volume and the high intensity focused ultrasound 
transducer array 12. Secondarily at risk are those tissues on 



ring to FIG. 8, at one step 82, diagnostic imaging is 
performed to identify and locate a treatment volume to be 
treated with HIFU energy. In doing so, diagnostic ultrasound 
pulses are transmitted from the imaging array portion 36 of 
the transducer array 12. The reflected ultrasound energy then 
is processed to view an image of the patient. An operator 
moves the transducer 12 and adjusts the focus or steering of 
the imaging array portion 36 to locate the desired treatment 
volume 16. At another step 84 a focused ultrasound beam is 



the other side of the treatment volume from the transducer 10 transmitted into the patient through a first aperture including 



array 12. 

Tne desired effect of treatment is heating of tissue within 
the treatment volume 16 resulting from absorption of the 
high intensity ultrasound. Some of the energy absorbed is 
received from the propagating ultrasound beam at the 
insonicating frequency. Additional energy is absorbed from 
ultrasound propagating at harmonics of the insonicating 
frequency. Such harmonics occur because of the nonlinear 
propagation of the ultrasound beam. This causes ultrasound 
energy at the insonicating frequency to be converted to 
ultrasound energy at harmonic multiples of the insonicating 
frequency. Higher frequency ultrasound is absorbed at a 
higher rate, proportional to frequency. Thus, there is an 
added and beneficial heating effect from the ultrasound 



20 



at least the imaging array portion 36 of the transducer array 
12. In one embodiment the beam is transmitted as a burst at 
a prescribed frequency for less than 5 microseconds in 
duration, resulting in a moderate intensity beam, (e.g., 1-20 
W/cm 2 SPTP, where SPTP means Spatial Peak Temporal 
Peak). The prescribed frequency is within the ultrasound 
frequency range and varies according to the application. The 
specific frequency is selected for the specific organ or tissue 
type being treated. 

A moderate intensity beam is used for the calibration step 
because a low intensity beam, such as used for the imaging 
step will not cause the desired harmonic generation in the 
focus zone that occurs using a moderate or high intensity 
ultrasound beam. Moderate intensity as used herein means 



beam harmonics. There is relatively greater absorption at 25 intensities at which the relationship between molecular 
these harmonics because of the greater absorption which velocity fluctuations and the pressure fluctuations within the 
occurs at the higher frequencies. tissue is nonlinear. Preferably the moderate intensity is 
To achieve maximum exposure at the treatment volume, chosen to generate harmonics only at the focus of the 
a maximum intensity beam, as limited by the aperture, is to ultrasound beam. For an imaging application a high con- 
be achieved. Such maximum intensity is to be achieved 30 ventional ultrasound pulse intensity is up to 10 W/cm 2 SPTP. 
without adversely impacting the tissue surrounding the At a decreased conventional imaging intensity of 100 
treatment volume. Dispersion of the ultrasound beam is one mW/cm 2 SPTP, the sinusoidal pressure fluctuation occurring 
source of difficulty. In particular refractive aberration due to in tissue is approximately 0.6 atmospheres (see waveform 
differences in the speed of ultrasound in the most superficial 29 of FIG. 1C). Thus, the highest instantaneous pressure in 
tissues nearest the ultrasound transducer cause the greatest 35 tissue is between 0.4 atm and 1.6 atm. A linear relationship 
dispersive effect on the focus at the treatment volume. This therefore exists between molecular velocity fluctuations and 
invention is directed toward a method which predicts refrac- pressure fluctuation. This is desirable for accurately imaging 
tive dispersion of the focus of the ultrasound on the treat- a treatment volume. At the threshold imaging intensity (i.e., 
menl volume, compensates for the predicted dispersion, and 300 mW/cm 2 SPTP), the sinusoidal pressure fluctuation 
tests the resulting focus to verify that the effects have been 40 occurring in tissue is approximately 1.0 atmospheres. Thus, 



the highest instantaneous pressure in tissue is between 0.0 
atm and 2.0 atm. Thus, a linear relationship still exists 
between molecular velocity fluctuations and pressure fluc- 
tuation. 

According to an aspect of this invention, the calibration is 
performed with a beam having an intensity outside such 
linear range in the focal zone even in the presence of 
superficial attenuation so as to benefit from the nonlinear 
response. The moderate intensity beam used initially is the 



properly compensated. As a result, exposure and harmful 
effect of the high frequency focused ultrasound on the 
surrounding tissue is minimized. 

Referring to FIGS. 6 and 7, the transducer array 12 
implements a HIFU treatment of a patient's organ 40. An 45 
acoustic coupling pad 13 is positioned between the patient's 
skin S and the transducer array 12. Ultrasound energy travels 
from the transducer elements 32, 34 through the patient's 
skin S and through various tissues 68 toward the organ 40. 
The ultrasound energy focuses at the treatment volume 16. 50 highest intensity available using the imaging array elements 
In an implementation where bone, air pockets or other 34. In later iterations during the calibration, the transducer 
ultrasound -opaque regions occur along the path between a elements 34 and 0 or more HIFU elements 32 define an 
transducer element and the treatment volume, such trans- aperture emitting ultrasound energy forming a beam of 
ducer elements are deactivated as shown in FIG. 7. moderate power intensity. An aberration will occur at such 
Specifically, FIG. 7 shows an ultrasound beam which travels 55 power level, so the delay lines 24 are adjusted to correct for 
through the intercostal regions of a patient's rib cage 72. such aberration (see step 92). 

According to an aspect of this invention, the aperture of the In an exemplary embodiment an attenuation derated spa- 
transducer array 12 is controlled so as not to transmit energy tial peak temporal peak intensity of approximately 10,000 
through the ultrasound-opaque regions. This is achieved by mW/cm 2 SPTP is the moderate power intensity achieved, 
monitoring the echoes of the transmitted ultrasound burst at 60 (i.e., about the highest ultrasound pulse intensities used for 



each active element. If the echo amplitude to that element, 
when gated to the target volume, is below a threshold value, 
then it is assumed that the path is blocked by an ultrasound- 
opaque region. Such elements then are deactivated. 

FIG. 8 shows a high level flow chart 80 for an embodi- 
ment of the method of this invention. FIGS. 9A-B show a 
more detailed flow chart for a similar embodiment. Refer- 



65 



diagnostic imaging and about 100 times less than HIFU 
(CW) pulse intensities). At such intensity the thermal and 
mechanical effects upon a patient are substantially less than 
for HIFU treatment. At the lower 300 mW/cm 2 intensity 
discussed above, the pressure variation is 1.0 atm. and the 
highest instantaneous pressure is between 0.0 atm and 2.0 
aim. Further increases in intensity cannot drive the lower 
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limit negative. As a result, at higher intensities (such as the and imaging transducer element 34 are adjusted to compen- 

moderate intensity beam of step 84) the pressure wave sate for phase change introduced by differing tissue along a 

becomes clipped, introducing a harmonic, including the path towards the treatment volume. Another cycle of the 

specific harmonic (e.g., the third harmonic) into the pressure time calibration then commences with the repcrformance of 

fluctuations. For such higher intensity pulses (e.g., 10,000 5 steps 84-92 until there is no phase error as determined at 

mW/cm 2 SPTP), echoes from the lower intensity regions step 90. 

will contain only the fundamental frequency, whereas ech- A time calibrated phase -inverted signal then is used for 
oes from the high intensity regions will contain the funda- transmitting the high intensity focussed ultrasound to the 
mental frequency and the specific harmonic frequency. The desired treatment volume. Thus, while phase inversion is 
reflected echoes from the treatment volume act as a point 10 applied, ultrasound energy is transmitted at step 94 from at 
source reflecting a nonlinear ultrasound burst (at the funda- least the HIFU array portion 38 to the treatment volume 16 
mental frequency and harmonics) to both the imaging ele- at a HIFU intensity for achieving desired ultrasonic medical 
ments 34 and the HIFU elements 32. Presence of the specific therapy of the treatment volume. Preferably a narrow band 
harmonic allows for measurement of the propagation lime transmit burst is used to achieve the desired harmonics at the 
from the treatment volume. In particular the nonlinear 35 treatment volume. The HIFU intensity is greater than the 
response is used as a signature to identify delay time for moderate intensity applied at step 84, (e.g., the HIFU 
echoes received at the various HIFU transducer elements 32 intensity is greater than 10 W/cm 2 SPTP). In one embodi- 
and imaging transducer elements 34. ment the HIFU intensity is at least 500 W/cra 2 SPCW. In an 
Accordingly, there is a threshold intensity level of the exemplary liver treatment application, a HIFU intensity of 
applied ultrasound energy, above which a relationship 20 1000 W/cm 2 SPCW was used (see waveform 21 of FIG. 
between the molecular velocity fluctuations and the pressure 1A). For an application as shown in FIG. 7, select elements 
fluctuations becomes nonlinear. Such nonlinearity is used to 32 of the HIFU array portion 38 overlaying a patient's ribs 
estimate the propagation time for ultrasound pulses trans- are deactivated while step 96 is performed to avoid bone 
mitted from respective imaging ultrasound elements 34 and heating, without degrading the effectiveness of other trans- 
received as echoes at HIFU transducer elements 32. In 25 ducer elements 32 emitting ultrasound toward the treatment 
particular, the echo at the specific harmonic of an ultrasound volume 16. According to preferred embodiments the pitch 
pulse of moderate intensity (e.g., 10 W/cm 2 SPTP) and less between transducer elements 32 is prescribed to be fine 
than HIFU intensity (1000 W/cm 2 SPCW), is used to enough to avoid phase aliasing and may be finer than the 
estimate the propagation path transit time for each HIFU spatial frequency of undulations of subcutaneous fat. At step 
transducer element 32 that will converge in a treatment 30 98, echoes of the HIFU beam are received at imaging 
volume 16 through nonhomogeneous tissue. The moderate elements 34 using harmonic bandpass filtering so as to 
intensity level 10 W/cm 2 SPTP is not enough to perform monitor the HIFU treatment. 

HIFU. Because the 10 W/cm 2 intensity is outside the linear FIG. 10 shows waveforms 100, 102, 104 which represent 

region between molecular velocity fluctuations and pressure the echo of transmitted ultrasound pulses at an exemplary 

fluctuations, third harmonic echoes are distributed in all 35 prescribed frequency of 0.95 MHz. Each waveform corre- 

directions from the treatment volume. These third harmonics sponds to a transmit signal at a different intensity. As 

echoes have a constant phase angle as the angle of scattering expected the echo drops off due to attenuation over time, 

is varied for each radial value (using a spherical coordinate wherein time corresponds to depth of reflection within the 

system centered on the target volume). The phase errors in phantom. FIG. 11 shows waveforms 200, 202, 204 which are 

this spherically expanding wavefront are of interest. 40 the third harmonic frequency (e.g., 2.85 MHz for the 0.95 

Specifically, the temporal delay in the specific harmonic MHz exemplary fundamental frequency) responses corre- 

echoes, such as the third harmonic echoes, observed by each sponding to the pulses of waveforms 100, 102, 104. FIG. 12 

given transducer element 32 serves as a measure of the shows the ratio of the corresponding waveforms in FIGS. 10 

propagation path transit time to transmit a pulse that will and 11 (i.e., the third harmonic magnitudes relative to the 

converge on the treatment volume 16. 45 fundamental frequency magnitudes). Referring to FIG. 12, a 

At step 86, echoes of the transmitted beam are received in first peak 132 at 30 microseconds echo delay corresponds to 
a second aperture including at least at the HIFU array a 23 mm depth. A second peak 134 at 47 microseconds echo 
portion 38, Due to the nonlinear relationship between delay corresponds to 36 mm depth. A third peak 136 at 62 
molecular velocity fluctuations and pressure fluctuations, the microseconds echo delay corresponds to 48 mm depth, 
echoes include ultrasound energy at the prescribed fre- 50 Peaks 132, 134 are locations in which near field peaks for the 
quency and at a specific harmonic of the prescribed fundamental frequency and third harmonic frequency over- 
frequency, such as the third harmonic. The specific harmonic lap. At the peak 136 the third harmonic echo from the 
echoes are of interest for determining phase advancement. highest intensity is much greater than the similar ratios for 

The specific harmonic echo resulting from the harmonic the other depths and from the lower transmit intensities. The 
conversion in the target volume due to the moderate inlen- 55 depth 48 mm for such peak 136 corresponds to the design 
sity beam is used to measure the relative propagation time focal depth (i.e., the treatment depth). Phase delay of one or 
for each HIFU transducer element 32. For embodiments in more elements 32 of the HWU array portion 38 and imaging 
which the elements 32 and 34 arc used during a HIFU elements 34 of imaging array portion 36 arc adjusted to alter 
treatment, the propagation time is determined for each the depth at which such maximum peak occurs (see step 92). 
element 34 also. At step 88, the quadrature pairs of the 60 The depth of the maximum peak is to be adjusted to the 
echoes are tested for phase error in the specific harmonic depth of the treatment volume to be treated, 
frequency, (e.g., third harmonic frequency). If there is no Referring to FIGS. 9A-B, the method is described in 
error, then HIFU treatment can commence (see step 90). If more detail. At one step 142 diagnostic imaging is per- 
there is an error, then at step 92 the delay lines 24 of the formed using the imaging array 36 to locate the treatment 
individual transducer elements 32, 34 are adjusted to offset 65 volume 16. At step 144 the transducer probe scanhead is 
dispersion effects of the nonhomogeneous tissue. More moved until the treatment volume 16 is located at a con- 
specifically, the phase of each HIFU transducer element 32 vergence point of HIFU array elements 32 at some depth D. 
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At slep 146 a first loop begins. At step 148 a variable W is times the quadrature component is less than the threshold 
set to 0. At step 150 an ultrasound pulse is transmitted from intensity, then transmit power is increased and the process 
the imaging array 36 to focus the transmit beam on the initiated again at the start of loop 1. Otherwise at step 198 
treatment volume 16. This is also the convergence point of phase angle error, 0, for such element 32/ is approximated 
the HIFU elements 32. At step 152 quadrature sampling is 5 as the quadrature component value divided by the in phase 
performed on each element 34/ of imaging array 36. At step component value or as the arc tangent of the quadrature 
154 four consecutive samples originating from the treatment component value divided by the in phase component value, 
volume 16 are processed for each element 34/. Such samples At step 200 the time delay for the element 32/ is adjusted 
are selected based upon the expected round trip time from a using the corresponding delay line 24 for the error in the 
respective element 34/ to the center of the treatment volume 10 time of flight using the fundamental frequency F. The time 
16 and back, including any time delay included at a corre- adjustment is the phase angle error for such element 34/ 
sponding delay line 24. At step 156 a second processing loop divided by the quotient 2 pi over 3 times the fundamental 
within the first processing loop begins. frequency. At step 202 a new time delay is derived for the 
At step 158 data for an element 34/ is selected for element 32/ by adding the adjustment to the prior value of 
processing. At step 160 the in phase component I is derived is the time delay for such element 32/. At step 204 another 
for the element 34/ by subtracting the third sample value element 32(i+l) is selected for processing. If all elements 32 
from the first sample value. At step 162 the quadrature have been processed then processing jumps to step 208, 
component Q is derived for the element 34/ by subtracting otherwise at step 206 processing goes back to the beginning 
the fourth sample value from the second sample value. At of the lo <>P 3 ( sle P 188 ) 10 process the next element 32(i+l). 
step 164, the results are tested to see if the threshold intensity 20 At ste P 208 the focussing is complete. At step 210 a HIFU 
has been exceeded. If the in phase component times the in hurst then is emitted to perform HIFU therapy on the 
phase component, plus the quadrature phase component treatment volume 16. 
times the quadrature component is less than the threshold Meritorious and Advantageous Effects 
intensity, then transmit power is increased and the process According to one advantage of the invention, HIFU 
initialed again at the start of loop 1. Otherwise at step 166 25 pulsesare transmitted through nonhomogeneous tissue with- 
phase angle error, E, is approximated as the quadrature out adversely heating tissue in the path between the ultra- 
component value divided by the in phase component value sound transducer and the treatment volume. Thus, treatment 
or as the arc tangent of the quadrature component value volumes can be treated to achieve hyperthermia, 
divided by the in phase component value. At step 168 the superthermia, ultrathermia or cavitation without adversely 
time delay for the element 34/ is adjusted using the corre- 30 impacting intervening tissue. 

sponding delay line 24 for the error in the time of flight using Although a preferred embodiment of the invention has 
the fundamental frequency F. The time adjustment is the been illustrated and described, various alternatives, modifi- 
phase angle error for such element 34/ divided by the cations and equivalents may be used. For example although 
quotient 2 pi over 3 times the fundamental frequency. At step a system in which each transducer element has a correspond- 
170 the adjustment is test. If the adjustment is less than zero 35 mg transmitter and receiver is described, in an alternative 
plus or minus some delta value, then the variable W is embodiment there is one transmitter which generates mul- 
incremented. At step 172 a new time delay is derived for the tiple drive signals for the respective transducer array ele- 
element 34/ by adding the adjustment to the prior value of ments. Although a preferred embodiment of the method has 
the time delay for such element 34/. At step 174 another been described for focussing a HIFU ultrasound burst upon 
element 34(i+l) is selected for processing. If all elements 34 40 a treatment volume, the method also may be applied for 
have been processed then processing jumps to step 178, focussing a diagnostic imaging ultrasound burst or beam at 
otherwise processing goes back to the beginning of the a target volume. Also, although in the preferred embodiment 
second loop (step 156) to process the next element 34(i+l). an imaging array and a HIFU array are shown and described, 
At step 178 the first loop's processing continues. At step the method also may be implemented using merely one 
180, W is tested to determine that is zero. If not 0 then 45 imaging transducer clement and/or only one HIFU trans- 
processing goes to the beginning of loop 1 (step 146). If duccr clement. Therefore, the foregoing description should 
W«0, then at step 182 an ultrasound pulse is transmitted not be taken as Urn iting the scope of the inventions which are 
from the imaging array 36 and one or more elements 32 of defined by the appended claims, 
the HIFU array 38. At step 184 quadrature sampling is What is claimed is: 

performed on each element 32/ of the HIFU array 38. At step 50 1 A method for applying ultrasound energy through a 

186 four consecutive samples originating from the treatment nonhomogeneous mass to an identified treatment volume to 

volume 16 are processed for each element 32/. Such samples provide ultrasonic therapy to the treatment volume, wherein 

are selected based upon the expected round trip time from a application of ultrasound energy causes molecular velocity 

respective element 32/ to the center of the treatment volume fluctuation and pressure fluctuation along a path from a 

16 and back to the element 32/, including any time delay, G, 55 source of ultrasound energy to the treatment volume, and 

included at a corresponding delay line 24. At step 188 wherein there is a threshold intensity level at the treatment 

another processing loop (i.e., loop 3) within the first pro- volume of the applied ultrasound energy, above which a 

ccssing loop begins. relationship between the molecular velocity fluctuations and 

At step 190 data for an element 32/ is selected for the pressure fluctuations becomes nonlinear in the treatment 

processing. At step 192 the in phase component I is derived 60 volume, the method comprising the steps of: 

for the element 32/ by subtracting the third sample value configuring a first transducer array to define a first 

from the first sample value. At step 194 the quadrature aperture, wherein the first transducer array is formed by 

component Q is derived for the element 32/ by subtracting a plurality of transducer elements; 

the fourth sample value from the second sample value. At while the first transducer array is configured to define the 

slep 196, the results are tested to see if the threshold intensity 65 first aperture, transmitting a pulse of ultrasound energy 

has been exceeded. If the in phase component limes the in from the first transducer array toward the treatment 

phase component, plus the quadrature phase component volume at a prescribed frequency and at a first peak 
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intensity greater than the threshold intensity level at the 
treatment volume; 
receiving echoes of the transmitted ultrasound energy at 
the first transducer array and at a second transducer 
array, the second transducer array having a plurality of 
transducer elements, wherein the echoes include ultra- 
sound energy at a specific harmonic of the prescribed 
frequency; 

for echoes received during an expected time window for 
receiving echoes from the treatment volume, compar- 
ing arrival time phase of the received echoes at the 
specific harmonic for each active element of the first 
transducer array and the second transducer array; 

adjusting phase delay for received echoes for one or more 
elements of the second transducer array to achieve 
phase inversion; and 

while phase inversion is applied to the second transducer 
array, transmitting ultrasound energy from the second 
transducer array to the treatment volume at a second 
intensity for achieving desired ultrasonic medical 
therapy of the treatment volume, wherein the second 
intensity at the treatment volume is greater than the first 
peak intensity at the treatment volume. 

2. The method of claim 1, wherein the steps of 
transmitting, receiving and comparing are repeated until 
each active element of the second transducer array is opti- 
mally aimed at the treatment volume. 

3. The method of claim 1, wherein the step of transmitting 
from the first ultrasound transducer array comprises emitting 
a pulse of less than 10 microseconds in duration. 

4. The method of claim 1, wherein the first peak intensity 
is at least 500 mW/cm 2 SPTP. 

5. The method of claim 1, wherein the first intensity is 
greater than 500 mW/cm 2 SPTP and the second intensity is 
at least 50 W/cm 2 SPTP. 

6. The method of claim 1 applied to a patient, wherein the 
treatment volume is in an organ, and further comprising 
prior to the step of transmitting from the first transducer 
array, the steps of: 

measuring amplitude of the received echoes; 

deactivating first elements of the first transducer array 
which receive echoes having an amplitude less than a 
threshold amplitude, wherein the threshold amplitude 
serves to identify ultrasound-obstructed pathways; and 

deactivating second elements of the second transducer 
array which receive echoes having an amplitude less 
than the threshold amplitude. 

7. A method for applying ultrasound energy through a 
nonhomogeneous path to an identified treatment volume to 
provide ultrasonic therapy to the treatment volume, wherein 
application of ultrasound energy causes molecular velocity 
fluctuation and pressure fluctuation along a path from a 
source of ultrasound energy to the treatment volume, and 
wherein there is a threshold intensity level of the applied 
ultrasound energy at the treatment volume, above which a 
relationship between the molecular velocity fluctuations and 
the pressure fluctuations becomes nonlinear, the method 
comprising the steps of: 

configuring a transducer array to define a first aperture, 
wherein the first transducer array is formed by a first 
plurality of imaging elements and a second plurality of 
treatment elements, the treatment elements capable of 
transmitting higher ultrasound power than the imaging 
elements, and wherein the first aperture includes at 
least multiple ones of the first plurality of imaging 
elements; 
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while the transducer array is configured to define the first 
aperture, transmitting a burst of ultrasound energy for 
a prescribed time from the transducer array toward the 
treatment volume at a prescribed frequency and at a 
first peak intensity greater than the threshold intensity 
at the treatment volume; 

configuring the transducer array to define a second first 
aperture, wherein the second aperture includes at least 
multiple ones of the first plurality of imaging elements 
and multiple ones of the second plurality of treatment 
elements; 

while the transducer array is configured to define the 
second aperture, receiving echoes of the transmitted 
ultrasound energy at the transducer array, wherein the 
echoes include ultrasound energy at harmonics of the 
prescribed frequency; 

for echoes received during an expected time window for 
receiving echoes from the treatment volume, compar- 
ing arrival time phase of the received echoes at the 
prescribed frequency and the specific harmonic of the 
prescribed frequency for each active element of the 
transducer array; 

adjusting phase delay for the received echoes for one or 
more of said elements of the transducer array to achieve 
phase inversion; 

configuring the transducer array to define a third aperture, 
wherein the third aperture includes at least multiple 
ones of the second plurality of treatment elements; and 

while phase inversion is applied to the transducer array, 
transmitting an ultrasound burst from the transducer 
array through the third aperture to the treatment 
volume, wherein the ultrasound burst is of a second 
peak intensity at the treatment volume for achieving 
desired ultrasonic medical therapy of the treatment 
volume, wherein the second peak intensity at the treat- 
ment volume is greater than the first peak intensity at 
the treatment volume. 

8. The method of claim 7, wherein the steps of configuring 
a first aperture, transmitting, configuring a second aperture, 
receiving and comparing are repealed until each active 
element of the second transducer array is optimally aimed at 
the treatment volume. 

9. The method of claim 7, wherein the step of transmitting 
from the ultrasound transducer array comprises emitting a 
pulse of less than 10 microseconds in duration. 

10. The method of claim 7, wherein the first peak intensity 
is at least 500 m W/cm 2 SPTP. 

11. The method of claim 7, wherein the first peak intensity 
is greater than 500 mW/cm 2 SPTP and the second peak 
intensity is at least 50 W/cm 2 SPTP. 

12. The method of claim 7, applied to a patient, wherein 
the treatment volume is an organ, and further comprising 
prior to the step of transmitting from the traasducer array, 
the steps of: 

measuring amplitude of the received echoes; 

deactivating elements of the transducer array which 
receive echoes having an amplitude less than a thresh- 
old amplitude, wherein the threshold amplitude serves 
to identify ultrasound -obstructed pathways. 

13. An ultrasound system for applying ultrasound energy 
through a nonhomogeneous path to a treatment volume to 
provide ultrasonic therapy to the treatment volume, wherein 
application of ultrasound energy causes molecular velocity 
fluctuation and pressure fluctuation along a path from a 
source of ultrasound energy to the treatment volume, and 
wherein there is a threshold intensity level of the applied 
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ultrasound energy at the treatment volume, above which a adjustable to introduce a select time delay respectively 

relationship between the molecular velocity fluctuations and for transmission and reception at a corresponding trans- 

the pressure fluctuations becomes nonlinear, the apparatus ducer element; 

comprising: means for defining a first aperture of the transducer array 
a first transducer array configured to define a first 5 through winch a burst of ultrasound energy is trans- 
aperture, the first transducer array having a plurality of mitted toward the treatment volume for a prescribed 
transducer elements which transmit ultrasound energy time at a prescribed frequency and at a first peak 
at a prescribed frequency and at a first peak intensity intensity greater than the threshold intensity at the 
greater than the threshold intensity at the treatment treatment volume, the first aperture including at least 
volume, the first aperture being formed by active ele- ]Q multiple ones of the first plurality of imaging elements; 
ments of the first transducer array; means for defining a second aperture of the transducer 
a second transducer array having a plurality of active array through which echoes of the transmitted ultra- 
transducer elements which receive echoes of the trans- sound energy are received, wherein the echoes include 
mitted ultrasound energy at a specific harmonic of the ultrasound energy at harmonics of the prescribed 
prescribed frequency; ^ frequency, the second aperture including at least mul- 
means for comparing phase of the echoes at the specific tiple ones of the first plurality of imaging elements and 
harmonic for each active element of the first transducer multiple ones of the second plurality of treatment 
array and the second transducer array by measuring elements; 

transit lime for the transmitted ultrasound energy to means for testing for time error of the echoes at the 

propagate from the first transducer array to the treat- ?Q prescribed frequency and the specific harmonic of the 

ment volume and echo back to the second transducer " prescribed frequency; 

array while the second transducer array is aimed at the means for reducing time error by adjusting respective 

treatment volume; and time delay at the plurality of time delay lines to achieve 

means for adjusting the phase of received echoes for one phase inversion and aim the transducer array at the 

or more elements of the second transducer array so that 2 s treatment volume; and 

each active element of the second transducer array is wherein while phase inversion is applied to the second 

aimed at the treatment volume, wherein propagation transducer array, ultrasound energy is transmitted 

transit time is applied as an advance time to each active through a third aperture of the transducer array to the 

element of the second transducer array to achieve phase treatment volume at a second peak intensity at the 

inversion; and 30 treatment volume for achieving desired ultrasonic 

wherein while phase inversion is applied to the second medical therapy of the treatment volume, wherein the 

transducer array, ultrasound energy is transmitted from second peak intensity is greater than the first peak 

the second transducer array to the treatment volume at intensity at the treatment volume, and wherein the third 

a second peak intensity for achieving desired ultrasonic aperture includes at least multiple ones of the second 

medical therapy of the treatment volume, wherein the 35 plurality of treatment elements, 

second peak intensity is greater than the first peak 18. The system of claim 17, further comprising means for 

intensity at the treatment volume. measuring propagation transit time for the ultrasound energy 

14. The system of claim 13, wherein a transmission of to propagate from the transducer array to the treatment 
ultrasound energy from the first ultrasound transducer array volume and echo back to the transducer array, and wherein 
comprises a pulse of less than 10 microseconds in duration. 40 adjustment of respective time delay is achieved by applying 

15. The system of claim 13, wherein the first peak the propagation transit times an advance time for each 
intensity is at least than 500 mW/cm 2 and the second peak element of the transducer array. 

intensity is at least 50 W/cm 2 SPTP. 19. The system of claim 17, wherein a transmission of 

16. The system of claim 13, further comprising a plurality ultrasound energy from the first ultrasound transducer array 
of time delay lines coupled to the transducer array, each of 45 comprises a pulse of less than 10 microseconds in duration, 
the time delay lines being adjustable to introduce a select 20. The system of claim 17, wherein the first intensity is 
time delay for the transmission at a corresponding trans- at least 500 mW/cm 2 SPTP and the second intensity is at 
ducer element of the transducer array, and wherein the phase least 50 W/cm 2 SPTP. 

adjusting means adjusts the time delay of at least one time 21. A method for applying ultrasound energy through a 

delay line to adjust the phase of the received echoes. 50 nonhomogeneous mass to an identified treatment volume to 

17. An ultrasound system for applying ultrasound energy provide ultrasonic therapy to the treatment volume, wherein 
through a nonhomogeneous path to a treatment volume to application of ultrasound energy causes molecular velocity 
provide ultrasonic therapy to the treatment volume, wherein fluctuation and pressure fluctuation along a path from a 
application of ultrasound energy causes molecular velocity source of ultrasound energy to the treatment volume, and 
fluctuation and pressure fluctuation along a path from a 55 wherein there is a threshold intensity level at the treatment 
source of ultrasound energy to the treatment volume, and volume of the applied ultrasound energy, above which a 
wherein there is a threshold intensity level of the applied relationship between the molecular velocity fluctuations and 
ultrasound energy at the treatment volume, above which a the pressure fluctuations becomes nonlinear in the treatment 
relationship between the molecular velocity fluctuations and volume, the method comprising the steps of: 

the pressure fluctuations becomes nonlinear at the treatment 60 transmitting a pulse of ultrasound energy from a first 

volume, the apparatus comprising: transducer toward the treatment volume at a prescribed 

a transducer array including a first plurality of imaging frequency and at a first peak intensity greater than the 

elements and a second plurality of treatment elements, threshold intensity level at the treatment volume; 

the treatment elements capable of transmitting higher receiving echoes of the transmitted ultrasound energy at 

ultrasound power than the imaging elements; 65 the first transducer and at a second transducer, wherein 

a plurality of time delay lines coupled to the transducer the echoes include ultrasound energy at a specific 

array, each one of the plurality time delay lines being harmonic of the prescribed frequency; 
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for echoes received during an expected time window for 
receiving echoes from the treatment volume, compar- 
ing arrival time phase of the received echoes at the 
specific harmonic for the first transducer and compar- 
ing arrival time phase of the received echoes at the 5 
specific harmonic for the second transducer; 

adjusting phase delay of received echoes for the second 
transducer to achieve phase inversion; and 

while phase inversion is applied to the second transducer, 
transmitting ultrasound energy from the second trans- 10 
ducer to the treatment volume at a second intensity for 
achieving desired ultrasonic medical therapy of the 
treatment volume, wherein the second intensity at the 
treatment volume is greater than the first peak intensity 
at the treatment volume. 15 

22. A method for focussing ultrasound energy through a 
nonhomogencous mass to a treatment volume, the method 
comprising the steps of: 

transmitting a pulse of ultrasound energy from a trans- 2Q 
ducer toward the treatment volume; 

receiving echoes of the transmitted ultrasound energy at 
the transducer, wherein the echoes include ultrasound 
energy at a specific harmonic of the prescribed fre- 
quency; 25 

for echoes received during an expected time window for 
receiving echoes from the treatment volume, compar- 
ing arrival time phase of the received echoes at the 
specific harmonic; 

adjusting phase delay of received echoes for the trans- 30 
ducer to achieve phase inversion; and 

while phase inversion is applied to the transducer, trans- 
mitting ultrasound energy from the transducer to the 
treatment volume. 

23. The method of claim 22, further comprising the step 35 
of receiving echoes of the ultrasound energy while phase 
inversion is applied to perform medical diagnostic ultra- 
sound imaging. 

24. 'Jne method of claim 22, wherein the transmitting 
ultrasound energy to the treatment volume while phase 40 
inversion is applied is for achieving ultrasonic medical 
therapy of the treatment volume in which application of the 
ultrasound energy causes molecular velocity fluctuation and 
pressure fluctuation along a path from the transducer to the 
treatment volume, and wherein there is a threshold intensity 



18 

level of the applied ultrasound energy at the treatment 
volume, above which a relationship between the molecular 
velocity fluctuations and the pressure fluctuations becomes 
nonlinear in the treatment volume. 

25. An ultrasound system for applying ultrasound energy 
through a nonhomogeneous path to a treatment volume to 
provide ultrasonic therapy to the treatment volume, wherein 
application of ultrasound energy causes molecular velocity 
fluctuation and pressure fluctuation along a path from a 
source of ultrasound energy to the treatment volume, and 
wherein there is a threshold intensity level of the applied 
ultrasound energy at the treatment volume, above which a 
relationship between the molecular velocity fluctuations and 
the pressure fluctuations becomes nonlinear, the apparatus 
comprising: 

a first transducer which transmits ultrasound energy at a 
prescribed frequency and at a first peak intensity 
greater than the threshold intensity at the treatment 
volume, the transducer and a second transducer receiv- 
ing echoes of the transmitted ultrasound energy at a 
specific harmonic of the prescribed frequency; 

means for comparing phase of the echoes at the specific 
harmonic for the first transducer by measuring transit 
time for the transmitted ultrasound energy to propagate 
from the first transducer to the treatment volume and 
echo back to the first transducer; 

means for comparing phase of the echoes at the specific 
harmonic for the second transducer by measuring tran- 
sit time for the transmitted ultrasound energy to propa- 
gate from the first transducer to the treatment volume 
and echo back to the second transducer; 

means for adjusting the phase of received echoes for the 
second transducer so that the transducer is aimed at the 
treatment volume, wherein propagation transit time is 
applied as an advance time to the second transducer to 
achieve phase inversion; and 

wherein while phase inversion is applied to the second 
transducer, ultrasound energy is transmitted from the 
second transducer to the treatment volume at a second 
peak intensity for achieving desired ultrasonic medical 
therapy of the treatment volume, wherein the second 
peak intensity is greater than the first peak intensity at 
the treatment volume. 

* * * * * 
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